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I. INTRODUCTION 

The crystalline phase breaks two continuous symme- 
tries: translational invariance and rotational invariancc. 
Dislocations are the topological defect associated with 
broken translational symmetry. Rotational distortions 
in crystals relax into wall structures formed from arrays 
of dislocations^ It would seem natural to describe the 
formation and evolution of these mesoscale dislocation 
structures using a continuum dislocation density theory. 
Individual dislocations are associated with long-range 
stress fields, and the dislocation evolution and struc- 
ture formation is strongly constrained by the need to 
screen these stresses. In the 1950s, a number of au- 
thorsi^^^ (building on ideas from differential geom- 
etry) developed such a continuum description involving 
the coarse-grained net topological charge of the disloca- 
tions in a region, organized into dislocation density ten- 
sor named for Nye. Since the long-range stress fields 
and the rotational distortions only depend upon the net 
dislocation density, it is natural for us to use this order 
parameter to describe the connections between rotations, 
dislocations, and stress. In this manuscript, we provide 
a systematic mathematical analysis of the relations be- 
tween rotational deformations, domain walls, and dislo- 
cation stress within the framework of Nye's continuum 
dislocation density tensor. 

In equilibrium, a large crystal with boundary condi- 
tions imposing a rotational deformation will relax by 
forming a few grain boundaries — sharp walls formed 
by dislocation arrays, separating perfect crystalline 
grains, with elastic stress confined locally near the grain 
boundaries.— In the continuum theory, the region of 
local elastic stress vanishes, and grain boundaries are 
described as (5-function dislocation densities with zero 



stress. We analyze stress-free dislocation walls in the 
continuum theory in section [PHI and connect them with 
Frank's formula^ for the dislocation content of grain 
boundaries in appendix [A"l 

At high temperatures, an initially disordered or micro- 
scopically deformed material will approach equilibrium 
through the formation and coarsening of polycrystals — 
the rotation gradients are confined to sharp grain bound- 
aries, separating local regions of different orientations. 
Again, the strain fields in a polycrystal are confined to 
small regions around the grain boundaries. In the contin- 
uum dislocation theory a polycrystal is thus a stress-free 
configuration of dislocations. In section [III we derive the 
most general solution for the dislocation density tensor 
with zero stress. In section IIII1 we show that this so- 
lution can be decomposed into a superposition of flat 
grain boundaries. In appendix [Bl and figure Owe explic- 
itly represent a (zero-stress) curved grain boundary as a 
continuum superposition of flat walls. In section HVl we 
also show that the most general stress-free state can be 
represented in terms of a local rotation field. Hence the 
stress-free states in the continuum theory can be inter- 
preted as polycrystals with arbitrarily small crystallites; 
the strain energy for grain boundaries is zero, so the crys- 
talline axes can vary arbitrarily in space. This would sug- 
gest that sharp, discrete wall formation is not implied by 
the energetics within the continuum theory; recent work 
elsewhere^ has shown that they may nonetheless form 
via shock formation in the natural dislocation dynamics. 

At low temperatures, where mass transport by diffu- 
sion of vacancies and interstitial diffusion is frozen out, 
external strain is relieved by (volume-preserving) dis- 
location glide. The motion, entanglement, and mul- 
tiplication of dislocations under this low-temperature 
plastic deformation leads to work hardening and 
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the development of a yield stress. These systems 
have long been modeled using continuum dislocation 
theories , 8 ' 9 ' 10 i 11 ' 12 i 13 ' 14 i 15 ' 16 i 17 sometimes simplifying the 
dislocation density into a scalar quantity, but some- 
times incorporating information not contained in the 
Nye tensor ('geometrically unnecessary' dislocation den- 
sities with canceling topological charge, yield stress laws, 
and separate dislocation densities for each slip sys- 
tem, ...). Continuum dislocation theories have been 
important not only in understanding large-scale dis- 
crete dislocation simulations) 1 ^ 18 ! 19 - 20 - 21 ^ 2 ! 23 ! 24 - 25 ! 26 - 27 
but also in understanding emergent dislocation avalanche 
phenomena^ 1 ^ 2 ^ 2 ^ 3 ^ 2 .^ 4 .^ At large deforma- 
tions (stage III plasticity), these tangles of dislocations 
begin to organize again into walls (here called cell walls) 
separating largely undeformed rcgion o 36 ' 37 ! 38 ' 39 ! 40 ' 41 with 
small misorientation angles between cells (again poten- 
tially explained by a shock formation in the climb-free 
dynamical evolution law 8 -^). In section [V] we provide 
the relation between the local rotation field for a gen- 
eral dislocation density tensor, which (for general stress- 
ful densities) is a formula providing a least-squares best 
approximation for the local orientation of the crystalline 
axes. 



II. THE GENERAL STRESS-FREE 
DISLOCATION DENSITY 

A dislocation is a crystallographic defect representing 
extra rows or columns of atoms and is characterized by 
two quantities; its line direction t, and its Burgers vector 
b. After a passage around a closed contour C that encir- 
cles a dislocation line, a displacement field u receives an 
increment b according to 

dui = j> pfj dxj = — bi (1) 

where (3 P describes the irreversible plastic deformation. 
The Nye dislocation density tensor p is defined by p = 
(t <g> b)<5(£) where £ is the two-dimensional radius vector 
measured from the axis of the dislocation in the plane 
locally perpendicular to t. When many dislocations la- 
beled by a are present, a coarse-graining description of a 
conglomerate of dislocations is preferred. In this picture, 
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with Gaussian weighting G(x — x') ~ 
(l/V2nL) 3 exp[-(x - x') 2 /(2L 2 )] over some dis- 
tance scale L large compared to the distance between 
dislocations and small compared to the dislocation 
structures being modeled. Since dislocation lines are 
topological and cannot end inside the crystal, p is 
divergence free: diPij = 0. 

A dislocation strains the crystal, and creates a long- 
range stress field. Peach and Kochler derived the rela- 
tionship for stress fields due to dislocations in an isotropic 



material^ 2 - In terms of the Nye dislocation tensor p, the 
stress can be written as the sum of two convolutions: 
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In Fourier space, the stress is given as a product: 

CT Q/3 (k) = Kap^Wp^Ck) , 
where the kernel 
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The problem of finding a family of dislocation config- 
urations with zero stress is equivalent to finding those 
densities p which are divergence free (ikiPij = 0) and are 
in the null space of K. Systematic investigation using 
Mathematica® shows that the solutions to the system of 
equations which incorporate both setting KijkmPkm = 
and i hi p^ = are 
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Since any given p has nine components and six con- 
straints, and since (e.g. for k = 0) K 1 has full rank, 
three basis tensors span the space of solutions. We solve 
for them and label them E x , E v and E z . 
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or simply 



Ejj — ik a 8ij + ikj8i a — iki£n m £j c 



(7) 
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Direct substitutions of the form of E a in place of p show 
that (j4| and the divergence free condition are simulta- 
neously satisfied for all values of a. It is convenient to 
include the imaginary number i into the expression for 
E a , because the Fourier transform of the gradient of a 
function is given by multiplying by ik. 

A general stress-free dislocation configuration there- 
fore can be written as a superposition of the three basis 
tensors E a : 



E«A a 



(9) 



The coefficients A z (k) form a valid vector field (i.e. A 
transforms like a vector). This vector will play a special 
role in determining the grain orientation inside each cell. 



Similarly, a twist boundary in the x-y plane can be 
generated by two sets of parallel dislocations oriented 
perpendicular to one another, one pointing in the x di- 
rection while another pointing in the y direction. It can 
be written simply as 



ptwist = _( 2vr )2 ^_ § (k x )6(k y )E* 
ik. 
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with the twist misorientation angle oj z = nb. The fact 
that one needs two perpendicular sets of parallel dislo- 
cations comes out naturally in this formulation. Because 
the number densities of the screw dislocations are the 
same in both directions, n here denotes the number den- 
sity in one of the two directions. 




III. DECOMPOSITIONS OF A STRESS FREE 
STATE INTO FLAT FRANK WALLS 



These three basis tensors can be used to describe grain 
boundaries. As an example, consider a tilt boundary in 
the x-y plane constructed from a set of parallel disloca- 
tion lines pointing along the x direction with the Burg- 
ers vector b pointing along the z direction. Let n be the 
number of dislocation lines per unit length along y. To 
make a plane in real space, we need two (5-functions in 
Fourier space. The boundary is then written 
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Notice that, for low-angle boundaries (small n), the tilt 
misorientation angle about the x axis is given by uj x = nb. 

We can write this tilt boundary in terms of our 
stress-free basis function E x . But why is the tilt- 
boundary stress free? Real grain boundaries have stresses 
from their constituent dislocations that cancel at long 
distances — they decay exponentially with distance over 
a length scale given by a typical distance between dislo- 
cations. Hence in the continuum limit where the disloca- 
tions become infinitely close together, the stress vanishes. 
Equivalently, the elastic energy of a boundary with low 
misorientation angle 9 goes as — 6 log 6*, which vanishes 
in the continuum limit b — > 0. Grain boundaries mediat- 
ing rigid rotations have vanishing stress in the mesoscalc 
continuum dislocation theory. 




FIG. 1: A general grain boundary whose normal is n posi- 
tioned at the distance A away from the origin separates two 
unstrained regions with a relative orientation defined by uj. 



A general boundary on the x-y plane is the sum of 
three types of boundary (a tilt along x, a tilt along y, 
and a twist along z): 
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where iv is the Rodrigues vector giving the angle of mis- 
orientation across the wall. The wall can be translated 
to a new position (z = A) by multiplying by e~ lk * A . The 
most general grain boundary with an arbitrary plane ori- 
entation can be obtained by then further rotating equa- 
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tion [12] by the rotation matrix 



R- 1 [n={0, ( j ) )} = [Ih{cj>) ■ Rp(-t 
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— sin(0) cos(0) 

sin(0) cos(0) 8111(6*) sin(0) cos(9) 



(13) 



to get 
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where il (9,(f>) defines a unit vector n normal to the 
plane of the boundary (see figure!]}, and k- A = kitnA = 
kiRi Z A = R~. fcjA. Equation[T4lis Frank's formula in the 
language of continuum dislocations. The connection with 
Frank's original formula is discussed in appendix [A"! 

To take this one step further, since it is possible to 
decompose any stress-free state into a linear combination 
of the tensor E a , it should also be possible to write a 
stress-free state as a superposition of flat cell walls. 



Theorem 1. Any stress-free state p SF can be written as 
a superposition of flat cell walls. Or more precisely, 
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where pf B is as previously defined, and 
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Proof. To get a general stress-free dislocation distribu- 
tion, one needs to integrate over three parameters de- 
noting the misorientation between the two grains, two 
angles defining each boundary, and the position of each 
grain component. 

To show this, we substitute the form of a[u>, Q, A] into 
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The integral over solid angle vanishes except along the 
line defined by the product of the two S- functions. Since 
our problem is isotropic, we may take this to be along the 
k x direction without loss of generality. Then the integral 
reduces to 

dcf> / <f (cos 9) 
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(18) 

where we use 6(g(x)) = J2 a ~ a )IW{ a )\i an( i the sum 
is taken over all a's with g(a) = and g'(a) ^ 042 The 

argument works for any k. Thus p? F (k) = E^A a (k) is 
shown. □ 

One must emphasize that this theorem does not ex- 
plain the prevalence of grain boundaries. Most stress- 
free states will be formed by continuous superpositions 
of walls. Indeed, even a curved grain boundary will de- 
mand such a continuous superposition (see appendix |B| . 



IV. STRESS-FREE STATES AND 
CONTINUOUS ROTATIONAL DEFORMATIONS 

In this section we show that the vector field A(x) intro- 
duced in the previous section is precisely the Rodrigues 
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vector field giving the rotation matrix that describes the 
local orientation of the crystalline axes at position x. 

What is A GB (x) associated with a grain boundary? 
Consider the form of A xy for a boundary lying in the x-y 
plane: 



A^(k) = (2tt) 



2 6( k x )S(k y ) 
ik z 



(19) 



The inverse Fourier transform of this expression involves 
an integral over a semi-circular contour in the upper com- 
plex plane, resulting in 

A^(x) = isign[z]u,. (20) 

In general, A GB is found after proper translation and 
rotation of the plane: 



A GB (x) = - sign[n ■ (x 



(21) 



The vector A(x) provides information about the local 
crystal orientation at the point x relative to a fixed global 
orientation. This is true in general: 

Theorem 2. The direction of A gives the axis of rota- 
tion of the local crystal orientation with respect to a fixed 
global coordinates by the amount provided by its magni- 
tude. 

In other words, the Rodrigues vector A(x) describes 
the local crystal orientations due to the presence of the 
stress-free dislocation density field p SF . 



Proof. First, note that 
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which, in real space, corresponds to 



pf? = <>,.\, - <),,(), „.\, 
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Now consider a rotation field Rij(A') = e e J' ifc k where 
A'(x) is the Rodrigues vector giving the local orienta- 
tion, and we wish to argue that A' can be used for A. 
Consider a small Burgers circuit C enclosing a region 5* 
with local orientation given by the field of A'(x). Inte- 
grating around the circuit C, the net closure failure — b 
due to the plastic distortion /3 P is given in terms of the 
local rotation -R(A') (see equation [1} : 

-bj= I Pl dx t = <f %(A'(x)) dxi (24) 
Jc Jc 

Applying Stokes' theorem to equation [24] and noting that 
the change in A' is small inside the small circuit C, we 
obtain 
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where we use the definition of the Nye tensor in the last 
equality. This expression holds regardless of the enclosed 
surface S, thus^i 
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□ 



Thus the stress- free distortions are precisely those gen- 
erated by rotation fields, and its dislocation density ten- 
sor field is given by our decomposition (equation l23[) with 
A equal to the Rodrigues vector for the local rotation. 



V. EXTRACTING THE LOCAL 
MISORIENTATION FROM THE NYE TENSOR 

The decomposition of jof? — A a Ef^ is somewhat differ- 
ent from the problem of breaking up a vector into pro- 
jections on various basis vectors. The main distinction 
lies in the fact that the three -E^-'s are not orthogonal 
to one another, so finding the components along them is 
not a simple dot product. We instead will minimize the 
square of the difference between the actual and the 

decomposition Ef- A". Let's define 
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Minimizing / will not only give the correct A Q for a 
stress-free p SF , it will also provide a natural definition 
for the local crystalline orientation of a general (stress- 
ful) dislocation density field. 

The minimization occurs when the derivative with re- 
spect to the component A 13 is zero: 
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It is possible to directly compute A in real space. From 

A = M i , Eij nn p mn 
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The expression of the Rodrigues vector A in real space, 
therefore by analogy to 1/k 2 factor in the Coulomb po- 
tential, is therefore 



A '<"> " 5 



(31) 



Equation [31] should be viewed as a natural definition of 
the local crystal axes, which could be invaluable for ex- 
tracting information about the misorientation angle dis- 
tribution, the wall positions, and hence the grain and cell 
size distributions^ 



VI. CONCLUSIONS 

In this manuscript, we explored the space of stress-free 
dislocation densities for an isotropic system. We showed, 
from first principles, that any stress-free state can be de- 
composed into a superposition of flat walls (grain bound- 
aries) and also can be written as a local rotational defor- 
mation field A(x). Finally, we provide a relationship be- 
tween this rotation field and the Nye dislocation density 
tensor, which in addition provides a formula for the best 
least-squares approximation for the rotation field for a 
stressful dislocation density. 

The analysis presented here forms the mathematical 
framework on which dynamical theories of continuum dis- 
location evolution are hung. It should offer basic tools 
for interpreting these simulations (identifying walls, mis- 
orientations, and rotational deformations during the evo- 
lution under polycrystallinc coarsening or plastic defor- 
mation), for theories based on the Nye tensor or more 
microscopic formulations. It should provide also a the- 
oretical basis for interpreting wall formation in contin- 
uum theories; minimizing stress provides a rationale for 
continua of walls, but not for discrete, individual grain 
boundaries or cell walls. 



APPENDIX A: FRANK'S FORMULA FOR A 
GENERAL GRAIN BOUNDARY: CONNECTION 
TO CONTINUUM THEORY 

Frank gave conditions on dislocation density for wall 
separating two perfect crystals mis-aligned by a rigid- 
body rotation^ Our analysis in section Hill explicitly gen- 
erated such walls, leading to a condition (equation IT3|) 



on the Nye dislocation density tensor. Here we relate 
Frank's original formulation with ours. For simplicity, 
we shall restrict ourselves to a small angle of misfit 9. 
For the treatment of large-angle boundaries, see Ref. [45| 
Let V be an arbitrary vector lying in the plane of a 
grain boundary, u> be an axis defining the relative rota- 
tion between the two grains separated by the boundary 
whose magnitude gives the net rotation angle 9, and b 
be the sum of the Burgers vectors of the dislocations cut 
by b, Frank's formula reads 

b = Vxu. (Al) 

(See Ref. |46| for the derivation^ and Ref. for the for- 
mula with an arbitrarily large angle 9.) 




FIG. 2: The orientation of the plain is defined by the vector 
normal n. The Rodrigues vector uj gives the axis of rotation 
and the angle of relative orientation between the two grains 
across the boundary. 

Using the Nye tensor, we can rephrase (|A1[) and then 
compare it with our statement of stress-free boundaries. 
Let's start off by defining a Burgers circuit C enclosing a 
surface S that intersects a grain boundary at two points 
x a and X;,. The net Burgers vector encompassed by the 
surface is b. Define V to be a vector lying in the bound- 
ary plane pointing from x a to x&, V = xj, — x a . We can 
represent this grain boundary by a constant matrix p° 
multiplied by a plane defined by 5(n ■ (x — A)), where n 
is a unit vector normal to the plane, and A is the per- 
pendicular vector pointing from the origin to the plane. 
(See figured) The integral of the Nye tensor p on the 
surface S gives the net Burgers vector b passing through 
that surface: 
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The (5-function serves to collapse the area integral into a 
line integral since the value is zero outside of the plane 
defined by x ■ n = 0: 



p% ; [nx V\idl 
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We can therefore relate the dislocation density to the 
rotation vector uj using Frank's formula (equation IA1|) : 
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With some relabeling, this becomes 

= (Pijh m + SijWm) ZimnVn - 
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Since V is an arbitrary vector in the plane of the grain 
boundary, we can write V as V = n x W for an arbitrary 
vector W. We can substitute n x W back into (IA5I). 
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This condition holds regardless of W. We can therefore 
safely ignore W in the equation. The condition now be- 
comes 
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The first term goes to zero because the first index of p?- 
designates the line component which always lies in the 
plane of the boundary. By definition, n is perpendicular 
to the plane, therefore, fiiP^j = 0. The condition for p° 
that makes a valid grain boundary is thus 

p°j = uj l n j - (u ■ n)<5y , 

or: 



(A8) 
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To see the connection between our formalism in ob- 
taining a general stress-free state, let us again rewrite the 
Fourier Transform of the general grain boundary p GB , 
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where all the variables are as defined previously. It is 
possible to perform the inverse transform of p GB to arrive 
at its real space representation. The two A-functions serve 
to define a plane in real space. The natural choice of 
coordinate is to make a rotational change of variables 
from (k x ,k y ,k z ) to (£ x , where & = R'/kj. In 



this coordinate, £ 2 is perpendicular to the plane of the 
boundary. The other two basis vectors lie in the plane of 
the boundary. 

The inverse transform can be written as 
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Note that since the new basis vectors are the rotation of 
the original set, its Jacobian is one. The next step is to 
express Ef- in terms of the new basis: 
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Similarly, 
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Substituting these into (|A10j) gives 
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The rotation matrix R was so constructed that Rz = n, 
or Ri Z = hi. Therefore, 

p% B = [uihj - (u ■ n)<%] 6(n ■ (x - A)) , (A14) 

exactly the same as what we derived from Frank's for- 
mula. 



APPENDIX B: DECOMPOSING STRESS-FREE 
STATES INTO FLAT WALLS: TWO EXAMPLES 

Here we illustrate theorem [T] and equation [H] by de- 
composing two stress- free states into a sum a[iv, CI, A] 
of flat Frank walls. Let's start with one of the simplest 
examples which is a flat twist boundary. According to 
equation [Til the boundary, in Fourier space, can be writ- 
ten as 



^twist 



nb 

ik z 



5{k x )8{k y )E z . 



GU) 
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The form of a[u>, CI, A], according to equation [TBI in this 
case is 



^ twist 



IUJ Z 



(27T) 3 7T 3 /: 

k 



dk' 



/3 {—nb) S(k' sin 9 cos cj>)5(k' sin 9 sin </>) ik > A 
ik' cos 9 



nbuj z e l w l <5(sin0 cos </>)<$ (sir 
(27r) 3 7r 3 /2 cos 6> 



2tt(5(A). 
(Bl) 



The combination of (^-function implies that <f> — or 4>, 
thus 

I |2 
7)/)/.; p~ M 

^ = ~ (2n)W ^( C0S ^( sin ^( A )> (B2) 

implying that such a wall can be created by only one 
regular straight wall. 




FIG. 3: (Color online) A circular grain boundary can be 

decomposed into a series of flat walls whose density decays as 
1/A 3 away from the center of the cylindrical cell. 



A more complicated example is the case where one 
cuts out a cylindrical portion of radius R inside a crystal 
with the axis of symmetry pointing along z, rotates it, 
and pastes it back (figure [3]). The resulting boundary is 
a circular grain boundary which can be represented in 
Fourier space as 



Ji 



k*R 



S(k z ) 



-.E z 



(B3) 



where J7i [•] is the Besscl function of type 1. In this case, 



iuj. 
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l«*>l 



dk' k' 



t Ji j\k'sm9\R) 
\k' sine | 



6(k' cos 9) e 
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5{cos9) / dfc'|fc'|Ji(|fc'|i?)e l ' c 



ik' A 



l2_„2//2 



2Ru>,e 



(27r) 3 7r 3 / 2 (A 2 -i? 2 ) 3 / : 



<5(cos6»)6(A- R). 



(B4) 



This example emphasizes the important point that we 
mentioned earlier, that a stress-free dislocation config- 
uration may need to be decomposed into a continuous 
superposition of flat cell walls. In particular, here we rep- 
resent a cylindrical wall as an infinite sum of flat walls 
with whose amplitudes go down as 1/A 3 with distance 
A away from the center of the cylinder. 
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